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Abstract 

Long-distance contributions of operators of the effective non-leptonic weak Hamilto- 
nian containing charmed currents have been recently studied. Penguin-like contractions 
of these operators, denoted as charming penguins, have been shown to be relevant for 
several B decays into two pseudoscalar mesons. In particular, they are expected to 
give large enhancements to processes which would be otherwise Cabibbo suppressed. 
Their contributions easily lead to values of BR{B~^ — > K^tt'^) and BR{Bd K~^7t~) 
of about 1 X 10~^, as recently found by the CLEO collaboration. In this paper, we 
show that such large branching fractions cannot be obtained without charming pen- 
guins. This holds true irrespectively of the model used to compute the non-leptonic 
amplitudes and of the relevant parameters of the CKM mixing matrix. We use the 
experimental measurements of the B~^ K^tt^ and B^ K^7r~ decay rates to con- 
strain the charming-penguin amplitudes and to predict BRs for a large set of other 
two-body decay channels where their contributions are also important. These include 
several pseudoscalar- vector and vector- vector channels, such as S — > pK and B — > pK*, 
which have not been measured yet. 
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1 Introduction 



In a recent paper [T| it has been shown that penguin-like contractions of operators 
containing charmed quarks, denoted as charming penguins, are able to enhance the 
S+ —f K^TT^ and Bd — >■ K~^tt~ decay rates with respect to the values predicted with 
factorization. By assuming reasonable values for the charming-penguin contributions, 
the corresponding branching ratios are of the order of (1-2) x 10~^, larger than the 
expected BR(Bd — > tt+tt"). This observation is particularly interesting because, in 
absence of charming-penguin diagrams, the B~^ — > K^Tr~^ and B^ — > K~^tt~ rates turn 
out to be rather small either because there is a Cabibbo suppression or because the 
non-Cabibbo suppressed terms come from penguin operators which have rather small 
Wilson coefficients a. The recent CLEO measurements Q 

BR{Bd^K+7r') = (1.5l[Ji;°i ± 0.1) X 10^5 

BR{B-^ ^ K\+) = {2.3tl:ltoi ± 0.2) x 10"^ 

BR{B+ ri'K+) = (7.lt|f ± 0.9) x 10"^ (1) 

BR{Bd v'K^) = (5-3^2:2 ± 1-2) x 10"^ 

allow an estimate of charming-penguin amplitudes and call for more quantitative stud- 
ies. 

In this paper, by using the experimental information obtained from the measured 
decay channels, we determine the parameters of charming penguins and predict a large 
set of BRs which have not been measured yet. The main results of our study, which 
will be described in detail below, are summarized in table In several cases we 
find that the value of the BR is strongly enhanced by charming penguins. The most 
interesting channels are Bd p~K+, B+ p+K° and B^ ojK^*^^, for which the 
contribution of charming penguins is larger than the theoretical uncertainties and the 
corresponding BRs are close to the present experimental upper limits. This means 
that they will eventually be measured, and compared to our predictions, in the near 
future. Within larger theoretical uncertainties, the BRs of other channels, such as 
B ttK*, B^ —5- ujK^*^^ and B — s- (pK^*\ are also appreciably enhanced by charming 
penguins and close to the present limits. Large enhancements are also expected for 
B — > r]^'^K decays. In these cases, however, our predictions are rather poor because 
of the presence of contributions, related to the anomaly, which are very difficult to 
evaluate. For B^, tt^tt^, we typically obtain a BR of about (5-10) x 10~^. Values as 
large as (2-3) x 10~^ remain, however, an open possibility. Finally, we also predict the 
BRs of several B — > Dh {h = vr, /?, etc.) decay channels, for which only upper bounds 
exist, see table 

The plan of the paper is the following: in section ^, we recall some basic facts about 
charming penguins and describe the main features of the present analysis; in section ^ 
we discuss the theoretical uncertainties stemming from the parameters of the CKM 
mixing matrix, the choice of the Wilson coefficients of the weak Hamiltonian and the 
models used in the calculation of the relevant amplitudes; in section |^ the procedure 
used in the determination of the important hadronic parameters is illustrated; sections ^ 
and m contain a brief discussion of the results and the conclusions respectively. 

^This is true unless the corresponding matrix elements are much larger than those predicted with 
factorization. 
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For the definition of tlie various parameters used in this study and the notation, the 
reader should refer to the experimental data have been taken from refs. ||2|, |3|, ^. 

2 Charming penguins 

"Charming penguins" denote penguin-like contractions in matrix elements of operators 
containing charmed quarks, the contribution of which would vanish using factoriza- 
tion H]. Their effects are enhanced in decays where emission diagrams are Cabibbo 
suppressed with respect to penguin diagrams. Similar effects can also be obtained 
by a breaking of factorization in the matrix elements of the penguin operators, Qs- 
QiQ, or assuming large chromo-magnetic contributions. Indeed all penguin-like effects, 
including the chromo-magnetic ones, have the same quantum numbers as charming 
penguins. For this reason they cannot be disentangled experimentally. Only an ex- 
plicit non-perturbative theoretical calculation of the operator matrix elements, done 
consistently in a given renormalization scheme (and missing to date), can distinguish 
the different terms. On the other hand, our parametrization of charming penguins 
effectively accounts for the other penguin effects. 

In absence of a consistent calculation of the matrix elements, penguin effects can 
also be enhanced by choosing a procedure where large "effective" Wilson coefficients of 
the penguin operators are obtained and factorization is used to compute the hadronic 
amplitudes. Attempts to reproduce the experimental data along this line have been 
recently made, see for example refs. js], ^ (the results of ref. Q will be discussed more in 
detail in section Q). In these approaches, the physical mechanism for the enhancement 
is similar to the one of ref. |l| , but the procedure used in the calculations is completely 
different: penguin diagrams involving charmed quarks are computed perturbatively on 
external quark states, and their effects included in the effective Wilson coefficients of the 
penguin operators Qs-Qio- It can be shown, however, that perturbation theory cannot 
be used to evaluate these effects: the relevant kinematical range of momenta involved 
in these calculations ^ corresponds to a region where charm-quark threshold effects, 
which are largely responsible for the enhancement of the B — > Ktt rates, are important 
and the perturbative approach ought to fail 0. For these reason, following ref. [|[, 
we parametrize charming-penguin effects in terms of the non-perturbative parameters 
riL and 6l, rather than make any attempt to compute them using perturbation theory 
and factorization. 

We now illustrate the main features of our new analysis: 

i) in ref. Q, all the amplitudes were normalized to the disconnected emission dia- 
grams, without computing them in any specific model. This makes the results 
model independent but only allows predictions of ratios of partial rates. In the 
present work we compute the disconnected-emission diagrams in the factorization 
approximation, which allows us to predict absolute branching ratios. The values 
of the hadronic parameters ^ and 6^, relevant for connected-emission diagrams, 
are extracted from a fit to decay channels which are expected to be dominated 
by emission diagrams, namely B ■kD^*\ B pD^*\ B K^*'iJ/^ and 
B —J- DD^*^; we also allow for some breaking of factorization, by varying the 
disconnected-emission amplitudes by 15% with respect to their factorized value; 
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ii) we do not fit the Wilson coefficients of the relevant operators of the weak Hamil- 
tonian, in the QCD combinations ai and 02, see for example ref. [^. We rather 
compute them at leading (LO) and next-to-leading (NLO) order in perturbation 
theory, at several values of the renormalization scale fx. In this way we estimate 
the uncertainty given by the scale and renormalization prescription dependence; 

iii) the values of the charming-penguin parameters rji and 6l are constrained by 
fitting the first two decay channels in eq. (|l]). We do not attempt to fit the B — > 
rj'K decay rates since these channels receive contributions, due to the anomaly, 
which are difficult to estimate. For these channels we will present in the following 
a rough estimate of the branching ratio, obtained by including non-anomalous 
contributions only. 

iv) using the values of the parameters as determined in i) and iii), a large set of 
BRs for two-body decay channels are predicted, including pseudoscalar-vector 
and vector- vector final states which were not considered in ref. |jl|. The results 
are also compared with the available experimental bounds; 

v) a comparative analysis of two-body BRs with and without charming penguins is 
performed. We show that there is strong evidence of a large charming penguin 
contribution in the data. As the experimental errors will be reduced, and more 
channels will be measured, this evidence will eventually find more support; 

vi) given its relevance in the extraction of the CP violating angle a [^, particular 
attention is devoted to BR{Bd tttt) decay. We will give an average value of 
the BR{Bd —5- vr'^vr'^), corresponding to central values of our parameters, and a 
maximum value, obtained by varying annihilation, GIM-penguin and charming- 
penguin contributions Q to this channel, with the constraint that the predictions 
of the measured channels remain compatible with the data at the 2-a level. 

In our analysis, we are mostly concerned with large enhancements induced by 
charming penguins, rather than with effects of the order of 20-30%. Moreover our 
study is not focused in testing the factorization hypothesis. For these reasons, since we 
obtain without any special effort satisfactory results for those decay channels which are 
dominated by emission diagrams, we have neither tried to optimize the parameters in 
this sector nor to estimate an error on the theoretical predictions. For these channels, 
we will only present results obtained with different models and Wilson coefficients. 
The spread of the results is representative of the theoretical error. 

On the other hand, the present experimental information on charming-penguin 
dominated decays is not sufficient to fix precisely the relevant non-perturbative pa- 
rameters. This is not a problem specific to our approach: in all cases where charming 
penguins are important, there are many smallish contributions (due to annihilations, 
GIM-penguins etc. |j|]) which cannot be precisely estimated by any non-perturbative 
method. Thus we will present a band of expected values for those channels where the 
enhancement is very large. 

3 Uncertainties of the theoretical predictions 

In order to assess the relevance of the charming penguins and make sensible predictions, 
a central issue is checking the stability of the results. We have studied their dependence 
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on several experimental and theoretical parameters which enter the calculation of the 
decay amplitudes. In particular: 

1. most of the theoretical predictions for two-body decays of heavy mesons are based 
on factorized formulae [l^-[12|. Although we do not assume factorization, we 



normalize the amplitudes to the factorized values of the disconnected-emission 
diagrams The latter depend on the particular model used to evaluate the 
relevant form factors [13|-||T^. We have selected a set of representative models 



which are consistent with the scaling laws derived from the HQET, namely lattice 
QCD (LQCD), the quark model of ref. |]l6| (QM) and the most recent results 
from QCD sum rules ||l^ (QCDSR). We also consider the model used in ref. ||T^ 
(ABLOPR), because a detailed analysis of the ratio BR{Bd — > K~^t:^) / BR[Bd 
vr+vr^) was presented there. We analyse the stability of the values of the fitted 
parameters rjL and 5l) and of the predicted BRs for different choices of 

the model used for computing the form factors; 

2. we vary the value of CKM-parameter a and the CP violating phase 5 (ue"*'^ = 
p — irj in the standard Wolfenstein parametrization We choose the ranges 

a = 0.24-0.48 and —0.7 < cos 5 < +0.7, although it can be argued on the basis 



of several analyses [|19|, |20[] that cos 5 > 0; 

3. we check the stability of our results against other effects which are suppressed in 
the factorization hypothesis. The contribution of annihilation diagrams and GIM 
penguins, as defined in ref. Q, has been studied. 

Before presenting our results, it is appropriate to discuss the evaluation of the 
relevant form factors used in the calculation of the disconnected emission amplitudes. 

Heavy-heavy matrix elements are related to the Isgur-Wise (IW) function ^(w). 
HQET allows the calculation of ^(w = 1) corresponding to q"^ ~ Qmax- most of the 
cases considered here, the range of momenta which is used is instead around q'^ ~ 0. 
Since the experimental measurements support a linear dependence of ^ on g^, the 
relevant parameter to evaluate the form factors is the slope of the IW function (a 
small extrapolation in is indeed needed since data exist down to rather small value 
of the momentum transfer). We find that the of our fit to the non-leptonic BRs is 
a steep function of p^. The data can be fitted reasonably well only if = 0.55-0.75 



(with a preferred value of 0.65), which lies within the allowed experimental range |21|. 
For more detailed and upgraded discussions on factorization and its tests, see ^. 

For the heavy- light form factors, HQET predicts scaling laws with the heavy quark 
mass only for q'^ ~ q'^ax @l • The evaluation of the form factors at the values of q^ 
needed for the factorized non-leptonic amplitudes requires further theoretical inputs, 
which are provided by a variety of non-perturbative approaches. The state of the art 
for the set of models used here is summarized in fig. |^. We find that the various 
determinations of the heavy-to-light form factors by the QM of ref. QCDSR |17|, 



and LQCD [15| agree reasonably well in a broad range of values of g^. For this reason, 
the values of the BRs obtained using these model are very similar, see tables |l|, |^ and 
^. On the contrary, ABLOPR systematically gives values of the form factors 1.5-2.5 
times larger than the other models. As a consequence most of the BRs computed with 
ABLOPR for B going into two light mesons are much larger than in all other cases. 
This fact has several consequences that will be discussed in sections ^ and ^. 
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Figure 1: Heavy-light form factors vs (GeV^) from LQCD (solid), QM (dashed), QCDSR 
(dot-dashed) and ABLOPR (hold- dashed) . 
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4 Fixing ^, S^, rji and Sl 



We now start the discussion of the results of our analysis. We stress again that, given 
the large uncertainties present in the evaluation of the BRs, our attitude is to focus 
mainly on large effects induced by charming-penguins. Although we do not even try 
to evaluate the theoretical errors, still we vary the input parameters, in order to check 
the stability of our predictions. 

The decay channels which have been used to fix the most important hadronic pa- 
rameters are presented in tables |3| and |^. In each table we show a list of decay channels 
calculated with and without charming penguins [in square brackets] B. The parameters 
of the charming-penguin contributions, i.e. rjL and 5l have been chosen as explained 
further on in this section. For each decay, we consider the results obtained using dif- 
ferent models for the form factors, namely the QM of ref. QCDSR0, LQCD m 



and ABLOPR |12|. The spread of the results obtained in the different cases can be 
taken as an estimate of the error due to the model dependence. These results have 
been obtained for "central values" of the parameters: the Wilson coefficients are com- 



puted at the LO at a scale ^ = 5 GeV (these coefficients are taken from ref. to 
which the reader can refer for details); cos 5 = 0.38 and a = 0.36; the disconnected 
emission diagrams are computed using factorization; GIM-penguin and annihilation 
contributions are put to zero. In the following the results obtained with this choice of 
the parameters will be denoted as CV (Central Values). 

Table |^ contains measured decay channels that are dominated by emission diagrams 
(and for which the dominant contributions are of 0(A^)), where A is the Cabibbo angle 
in the Wolfenstein parametrization. They are used for fitting ^ and 5^. Reasonable 
values of the parameters are obtained in all cases, with fairly good values. The 
channels dominating the fit, namely those with the smaller relative experimental error 
(the BRs of which are given above the horizontal line in the table), are essentially 
independent of the charming-penguin contributions and of the CKM parameters cos 5 
and a. For this reason, both the the fitted parameters have negligible depen- 

dence on these quantities. The experimental values of the BRs, as well as the values 
of the «^iid the results of the fit, for each model, are also given. For the different 
models, different determinations of ^ and 5^ are obtained. The larger ABLOPR form 
factors prefer an appreciably smaller value of ^ and have some problem in reproducing 
the measured BR{B K*J/^). 

In table ^ for completeness, we also present (for QCDSR only) a band of results 
obtained in the following way: we vary at the LO the renormalization scale from /i = 5 
io jJL = 2 GeV; we take the Wilson coefficients computed at the NLO in the NDR 
scheme (but only at /i = 5 GeV); we change the value of the disconnected-emission 
diagrams by ±15% with respect to their factorized value; we allow the value of the 
GIM-penguin and of the annihilation diagrams to vary within the intervals defined in 
ref. H and take -0.7 < cos (5 < +0.7 and 0.24 < o" < 0.48. The spread of values given 
in the table is representative of the theoretical uncertainty due to the parameters listed 
above. The band of results obtained by varying the parameters as explained above is 
denoted as the Band Value (BV). 

With the exception of ABLOPR, very similar results are obtained with all the 



^Only in those cases where charming-penguin contractions contribute. 
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models considered in our study. For this reason, in table ^ we only discuss the results 
obtained using QCDSR, as a representative case, and ABLOPR. In this table, we give 
a detailed "map" of the results obtained for BR{B+ K^tt+) and BR{Bd K+tt'), 
which have been used to fit the charming-penguin parameters ijl and 6l- We also 
present, for the same values of the parameters, BR{Bd tt+tt"). The values in square 
brackets are computed without charming penguins. To illustrate the effects of a large 
coefficient for the most important penguin operator Qq, we also give the results with the 
coefficients computed at the NLO in NDR: in this regularization scheme the relevant 
penguin coefficients turn out to be larger at NLO than at LO. The values of ^ and 6^ 
are taken from the fit to the data of table ^, which only involve heavy-heavy and heavy- 
light meson final states. For light-light meson final states, it would have been possible, 
of course, to choose values of ^ and 6^ different from those fitted on the heavy-heavy 
and heavy-light decay channels. Different choices of ^ for pseudoscalar-pseudoscalar, 
vector-pseudoscalar and vector-vector final states are also possible in principle. The 
inflation of parameters would have only complicated the analysis without changing the 
basic conclusions. 

When charming penguins are included, QCDSR (as well as with LQCD or the QM), 
give good results for BR{B^ — > K^tt~^) and BR{Bd K^tt^) and at the same time 
respect the experimental upper bound on BR{Bd tt+vt^). This is true almost ir- 
respectively of the choice of the Wilson coefficients and of cos 6 i and of a: only for 
large values of a, and very small values of cos 5, there is a problem with the upper 
limit on BR[Bii 7r"^7r~). There is a general tendency to have BR{B^ —>■ K^7:~^) 
close to BR{B(i — > i^+vr^) instead of larger, as found experimentally. Given the ex- 
perimental errors (and the uncertainties coming from GIM-penguins and annihilation 
diagrams) we believe that this difference is at present not significant. As far as the 
last of the BRs in eq. (2) is concerned, with charming penguins we typically obtain 
BR(B+ 7]'K+) ~ (2.5-3.0) x I0~^ slightly less than one-half of the central ex- 
perimental value (without charming penguins BR{B^ rj'K^) ~ (0.2-1.5) x 10~^). 
Since the difference probably comes from anomalous contributions, which are not di- 
rectly related to charming penguins and are very difficult to evaluate, this discrepancy 
does not affects our main conclusions. 

With charming penguins, ABLOPR may find even better agreement with the data, 
in that it easily reproduces the pattern BR{B~^ — > K'^tt^) ~ 1.5 BR{B(i K~^7r~) 
in a wide range of cos 5 and a. However, it violates the bound on BR{Bd tt^tt^) 
unless special values of the parameters are assumed. If we demand good agreement for 
all the three BRs, only a low value of a is acceptable. In this case, with a larger value 
of the coefficient of Qq, as obtained at the NLO, ABLOPR finds agreement with the 
data even without charming penguins. 

With the exception of ABLOPR, which we discuss separately below, we conclude 
that charming-penguin contributions are essential for reproducing the data. For all the 
other models that we have studied (QM, LQCD and QCDSR), theoretical estimates 
obtained using factorization without charming penguins predict too small BRs. In 
particular BR{B^ K^n^) is always a factor ~ 4-10 smaller than its measured 
value. Although the experimental errors are still rather large, we believe that these 
data already show evidence of substantial charming-penguin contributions to -B — > Kn 

•^Notice that with or without charming penguins BR{B^ K^tt^) is practically independent of cos (5. 
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decays. 

ABLOPR deserves a separate discussion. By using some specific set of Wilson 
coefficients (corresponding to a large value of the coefficient of the penguin operator 
Qq), this model may predict correctly the experimental values for BR{Bd K~^it~) 
and BR{B~^ — s- K^7r~^) even without charming-penguin contributions. This is the 
case, for example, if one computes the coefhcients at the NLO, or at a low value 
of the renormalization scale /i. The reason is that the values of the form factors in 
ABLOPR are quite atypical, see fig. |l]. The requirement that the experimental values 
of BR{Bd — > K~^'K~) and BR{B~^ K^7r~^) are reproduced without violating other 
experimental bounds, in particular BR{B(i 7r"'"7r~), demands, however, a special 
choice of the Wilson coefficients of the operators of the effective Hamiltonian and of 
the values of a, which must be small, or cos 5, which must be negative (see table |^. 
Other problems afflict this model, particularly in the vector sector. For example, many 
experimental upper bounds, see table ||, are violated by using ABLOPR. This model 
has also problems with semileptonic decays because it gives a very large semileptonic 
BR{B p) (and a relatively large BR{B ^ vr)): one finds BR{B p)/BR{B 
vr) = 2.98, to be compared to the experimental number 1.4]'^q4 it 0.3 it 0.4 |^] and, as 
far as Vub is concerned, Vuh/Vch = 0.047 from B ^ p and Vuh/Vcb = 0.070 from S — > vr. 

At present, we cannot completely exclude that, by using ABLOPR, it is possible to 
find a very special set of parameters satisfying the following requirement: large values 
of BR{Bd — > K~^7r~) and BR{B~^ — > K^ir'^), without invoking charming-penguin 
effects, and respect of the experimental bounds, notably those on BR{Bd — vr+vr^) 
and on BRs of channels containing vector mesons in the final states. Considering also 
the problems of ABLOPR with semileptonic decays, and the fact that it is so different 
from all the other models, we find this interpretation of the data rather remote. 

We now discuss the results of ref. Q where the B Kit, B rj'K and B — > vrvr 
branching ratios have also been studied. The authors of this paper conclude that 
is possible to describe the data by using a specific set of effective coefficients, which 
include charm-quark penguin effects computed in perturbation theory, and factorized 
hadronic amplitudes. 

On the theoretical validity of the approach followed in ref. [Q, we have already 
commented in section |2[ Concerning the phenomenological analysis, our main obser- 
vation is that, whenever (any kind of) penguins are important, the values of the BRs 
are strongly correlated, as shown in table ^. As an example, a very small value of ^, as 
required to reproduce BR{B^ ri'K~^), leads to a violation of the experimental bound 
on BR{B(i —5- vr^vr"). For this reason, in this table, we have decided to present the 
BRs for several choices of the main parameters separately. From the analysis of ref. , 
one may conclude that, for any given decay channel, there are points in the parameter 
space where agreement between theory and data can be achieved. However this is no 
longer the case when the correlations among different predictions, induced by using 
the same set of parameters for all the light-light channels, are taken into account. In 
particular, on the basis of our experience, we think that it is very difficult to reproduce 
the experimental BR{B~^ — > r]'K^) without violating the bound on BR{B(i — > -k'^tt"). 
Clearly, extra-contributions coming from the anomaly may change this conclusion. 
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5 Predictions for yet unmeasured BRs 



Having determined the hadronic parameters from the measured BRs, we use them to 
predict those which have not been measured yet. These include many channels where 
the charming-penguin contributions are Cabibbo-enhanced with respect to emission di- 
agrams, together with some Cabibbo-suppressed channels (of 0{X'^)), which are phys- 
ically interesting, such as B ^ vrvr and B vrr/'. The results are collected in table 0. 
As expected, charming-penguin-dominated channels are strongly enhanced: typically 
by a factor of 3-5; in some case, the enhancement can be as large as one order of 
magnitude. B pK^*\ B ujK^*^ and B rj^'^K^*^ are the channels most sensitive 
to charming-penguin effects. In some case, such as -B ^ r]K and B^ — > p^K^, the 
BRs remain small, at the level of 10~®, in spite of the large enhancement, and may be 
difficult to measure. 

By looking to the column QCDSR-BV, we observe that the typical uncertainty is 
about a factor of 2, although, for some channel, it can reach even factors of 10 or more. 
For this reason, for the Cabibbo-suppressed channels (corresponding to the BRs below 
the horizontal line) , the errors are so large that predictions with and without charming 
penguins largely overlap. This (helas!) occurs also in some charming-penguin-enhanced 
channels (the BRs of which are given above the horizontal line). All predictions in- 
volving rj' and rj suffer from the poor theoretical control that we have on anomalous 
contributions, which are never included in our calculations. 

The best candidates to test charming penguins are B — > pK decays (in the different 
charge combinations) and B^ ujK^*'^^. For these channels, the enhancement is larger 
than the estimated uncertainties, and the expected BRs are not far from the present 
experimental limits, so that one can expect that they will soon be measured. Within 
larger uncertainties, other channels are also appreciably enhanced by charming pen- 
guins, e.g. B ttK*, B+ LoK^*^+ and B (t)K^*\ With QCDSR, LQCD and the 
QM the experimental bounds are generally respected; in many cases, using ABLOPR, 
they are instead violated, often by a large factor, even in absence of charming-penguin 
contributions. 

The Cabibbo-suppressed B ^ tttt decay plays a fundamental role in the extraction 
of unitary-triangle angle a. In ref. Q, it has been shown that charming-penguin 
contributions in B^ 'k~^-k~ decays do not allow the extraction of this parameter using 
only the asymmetry measured in this channel. Still, it may be possible to extract a 
by using the isospin analysis of ref. |23]. This method requires BR{Bii ir^ir^) to 
be large enough. In ref. Q, by normalizing the rates to some conventional value for 
BR{Bii — > TT^vr^), it was shown that values of this BR as large as (2-3) x 10~^ were 
possible. In the present case, we typically obtain, within large uncertainties, an average 
value of (5-10) x 10^^. We then include other contributions, such as annihilation 
diagrams and GIM penguins [|^, and try to maximize the BR^B^ ir^ir^) as a function 
of the hadronic parameters, allowing x^-values twice larger than those found in the 
fit. In this case, by stretching all the parameters, we find that a maximum value 
BR{Bd vr^TT*^) ~ (2-3) x 10~® is still allowed. Therefore the extraction of a from the 
Bd 7r"'"7r~ CP-asymmetry may be a very difficult, hopefully not impossible, task. 

Finally, in table |^, our predictions of other emission-dominated channels, for which 
only upper bound have been measured, are also presented. We see that all the models 
satisfy the bounds without any problem. 
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6 Conclusions 



In this paper, by using the present experimental measurements and constraints, we 
have determined some non-perturbative hadronic parameters necessary to predict the 
BRs for two-body i?-meson decays. Using this information, we have predicted a large 
set of yet unmeasured BRs and discussed the corresponding uncertainties. The latter 
have been estimated by using different models, different Wilson coefficients and by 
varying the CKM parameters within the ranges allowed by experiments. 

We found that, in several cases, notably B pK, charming-penguin effects are 
larger than the theoretical error and the expected value of the corresponding BRs 
are close to the experimental upper bound. Therefore these channels are the ideal 
playground to test non-factorizable charming-penguin effects. 
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Channel 


QCDSR-CV 


QCDSR-BV 


Lattice-CV 


QM-CV 


ABLOPR-CV 


Experiment 




BR X 10^ 


BR X 10^ 


BR X 10^ 


BR X 10^ 


BR X 10^ 


BR X 10^ 


Bd ri'K" 


2.45 


[0.19] 


1.97-3.55 [0.12-0.67] 


2.49 


[0.10] 


2.39 [0.19] 


4.48 


[0.75] 




Bd r]'K*° 


0.75 


[0.01] 


0.30-2.41 [0.00-0.10] 


1.41 


[0.00] 


0.47 [0.01] 


2.56 


[1.06] 


< 9.9 


Bd r]K° 


0.07 


[0.01] 


0.01-0.26 [0.00-0.02] 


0.07 


0.00] 


0.06 [0.01] 


0.10 


0.02] 




Bd r]K*^ 


0.04 


[0.02] 


0.01-0.22 [0.01-0.06] 


0.05 


[0.01] 


0.03 [0.02] 


0.87 


[0.83] 


< 3.3 


Bd ujK^'' 


0.73 


[0.04] 


0.13-1.70 [0.03-0.38] 


0.80 


[0.02] 


0.70 [0.06] 


0.64 


[0.02] 




Bd ojK*° 


1.85 


[0.11] 


1.01-3.26 [0.08-0.58] 


2.76 


[0.12] 


1.43 [0.08] 


6.26 


[0.41] 


< 3.8 


B^ (fiK'^ 


1.48 


[0.11] 


0.68-2.29 [0.10-0.75] 


1.77 


[0.08] 


0.92 [0.12] 


2.59 


[0.92] 


< 4.2 


B^ (fiK'^ 


3.27 


[0.09] 


1.90-6.73 [0.07-0.62] 


5.08 


[0.11] 


1.65 [0.05] 


10.3 


[0.29] 


< 2.2 


B^ ttOrO 


0.65 


[0.04] 


0.49-0.95 [0.01-0.26] 


0.68 


[0.03] 


0.64 [0.03] 


0.94 


[0.13] 


< 4.1 


B^ -kOK-^ 


0.40 


[0.03] 


0.29-0.60 [0.01-0.13] 


0.61 


[0.02] 


0.26 [0.03] 


0.66 


[0.09] 


< 2.0 


Bd TT-K* + 


3.13 


[0.17] 


2.05-5.07 [0.05-1.01] 


3.47 


[0.11] 


3.16 [0.16] 


2.84 


[0.44] 


< 6.7 


Bd -* p-K+ 


0.63 


[0.02] 


0.21-0.67 [0.00-0.11] 


1.19 


[0.02] 


0.53 [0.02] 


1.89 


[0.09] 


< 3.3 


Bd -» p~K*+ 


2.21 


[0.12] 


1.44-3.57 [0.03-0.72] 


4.15 


[0.13] 


1.75 [0.09] 


7.57 


[1.18] 




Bd p°K° 


0.44 


[0.03] 


0.07-1.63 [0.00-0.09] 


0.81 


[0.01] 


0.35 [0.03] 


1.33 


[0.43] 


< 3.0 




1.01 


[0.01] 


0.50-2.45 [0.01-0.13] 


1.86 


[0.02] 


0.60 [0.01] 


4.39 


[0.16] 


< 46 


B+ ri'K+ 


2.31 


[0.20] 


2.05-3.14 [0.11-0.92] 


2.36 


[0.10] 


2.27 [0.20] 


3.80 


[0.70] 




B+ rj'K*+ 


4.14 


[0.04] 


1.45-11.22 [0.00-0.57] 


4.82 


[0.02] 


4.18 [0.03] 


85.0 


[1.59] 


< 29 


B+ riK+ 


0.12 


[0.02] 


0.02-0.46 [0.00-0.08] 


0.13 


[0.01] 


0.11 [0.01] 


0.25 


[0.05] 




B+ — > r]K*+ 


0.27 


[0.04] 


0.04-1.19 [0.01-0.22] 


0.29 


[0.03] 


0.25 [0.04] 


3.13 


[1.09] 


< 24 


B+ — » u}K+ 


0.60 


[0.02] 


0.06-1.39 [0.00-0.66] 


0.65 


[0.04] 


0.60 [0.03] 


0.39 


[0.12] 




B+ uiK*+ 


1.50 


[0.08] 


0.84-2.48 [0.03-1.12] 


2.24 


[0.16] 


1.20 [0.06] 


3.64 


[0.59] 


< 11 


B+ 4>K+ 


1.54 


[0.12] 


0.70-2.38 [0.10-0.78] 


1.84 


[0.09] 


0.96 [0.13] 


2.69 


[0.96] 


< 0.53 


B+ -» <I>K*+ 


3.35 


[0.09] 


1.96-6.90 [0.08-0.64] 


5.21 


[0.11] 


1.70 [0.06] 


10.6 


[0.30] 


< 4.1 


B+ n°K+ 


0.94 


[0.18] 


0.67-1.03 [0.09-0.74] 


0.86 


[0.10] 


0.95 [0.18] 


0.94 


[0.45] 


< 1.6 


B+ n'>K*+ 


1.70 


[0.14] 


1.15-2.42 [0.06-0.71] 


1.81 


[0.09] 


1.71 [0.13] 


1.55 


[0.35] 


< 8.0 


B+ ^ vr+X-o 


0.92 


[0.13] 


0.73-1.15 [0.08-0.37] 


1.29 


[0.07] 


0.62 [0.12] 


1.54 


[0.40] 


< 3.9 


B+ ^ p+K" 


2.67 


[0.00] 


0.88-7.28 [0.00-0.36] 


3.05 


[0.00] 


2.62 [0.00] 


27.3 


[0.00] 


< 6.4 


B+ ^ p+K"" 


2.21 


[0.09] 


1.42-4.01 [0.06-0.38] 


3.82 


[0.09] 


1.34 [0.07] 


9.50 


[1.06] 


- 


B+ ^ p"K+ 


0.39 


[0.06] 


0.18-0.61 [0.00-0.21] 


0.61 


[0.03] 


0.35 [0.06] 


0.81 


[0.53] 


< 1.4 


B+ ^ p"K*+ 


1.26 


[0.15] 


0.81-1.62 [0.06-0.63] 


2.18 


[0.14] 


1.00 [0.11] 


4.31 


[1.18] 


< 90 


B(i — > 7V~^7r~ 


0.99 


[0.63] 


0.39-2.63 [0.15-1.55] 


0.71 


[0.39] 


0.95 [0.60] 


2.37 


[1 711 


< 1.5 


Bd n°Tr° 




0.04 


0.01-0.12 [0.01-0.11] 


0.004 [0.002] 


0.06 [0.05] 




0.01 


< 0.93 


Bd p°-K° 


0.08 


[0.11] 


0.00-0.21 [0.02-0.21] 


0.01 


[0.07] 


0.12 [0.13] 


0.00 


[0.03] 


< 1.8 


Bri — > o+7r~ 


2.67 


1 71 


1.02-7.51 [0.45-4.15] 


1.93 


[1.07] 


2.55 [1.63] 


6.29 


[4.52] 


< 8.8 


Bd — ^ '^'^ P~ 


0.59 


[0.38] 


0.22-1.75 [0.11-0.91] 


0.66 


[0.37] 


0.46 [0.30] 


1.92 


[1.38] 




B . ^ „0 


0.16 


[0.18] 


0.00-0.58 [0.04-0.46] 


0.02 


[0.18] 


0.21 [0.19] 


0.10 


[0.19] 


< 28 


Bri — > 0~ 


1.81 


[1.16] 


0.69-5.10 [0.31-2.82] 


2.19 


[1.22] 


1.37 [0.88] 


16.6 


[11.9] 


< 220 


Bd — > wtt" 


0.06 


[0.01] 


0.03-0.42 [0.00-0.29] 


0.08 [0.005] 


0.06 [0.02] 


0.06 


[0.00] 




Bd — > i^p" 


0.10 


[0.01] 


0.05-1.22 [0.01-0.73] 


0.14 


[0.01] 


0.08 [0.01] 


0.36 


[0.01] 


< 3.4 


Bd -+ 4>TT° 


0.00 


[0.00] 


0.00-0.01 [0.00-0.01] 


0.00 


[0.00] 


0.00 [0.00] 


0.00 


[0.00] 


< 0.65 


Bd r]'-^° 


0.48 


[0.36] 


0.04-2.14 [0.01-2.91] 


0.21 


[0.22] 


0.44 [0.33] 


1.05 


[0.98] 


< 2.2 


B+ tt+tt" 




0.53 


0.16-1.18 




0.36 


0.51 




1.07 


< 2.0 


B+ p+tt" 




1.12 


0.30-2.62 




0.80 


1.05 




2.63 


< 7.7 


B+ TT+p" 




0.52 


0.19-1.08 




0.47 


0.47 




1.01 


< 5.8 


B+ p+pO 




0.94 


0.23-2.11 




1.08 


0.73 




7.44 


< 100 


B+ — > ui7r+ 


0.84 


[0.50] 


0.38-2.35 [0.08-1.64] 


0.94 


[0.45] 


0.70 [0.45] 


1.69 


[0.99] 




B+ u)p+ 


1.81 


[0.94] 


0.80-5.57 [0.19-3.40] 


2.44 


[1.04] 


1.36 [0.72] 


12.9 


[7.23] 




B+ <f>1T+ 


0.00 


[0.00] 


0.00-0.02 [0.00-0.02] 


0.00 


[0.00] 


0.00 [0.00] 


0.00 


[0.00] 


< 0.56 


B+ TT+V 


1.09 


[0.76] 


0.01-3.33 [0.13-3.36] 


0.87 


[0.58] 


0.97 [0.65] 


3.09 


[2.53] 


< 4.5 


Bd ^ ttO J/^- 




2.27 


2.06-2.36 




1.82 


2.70 




0.41 


< 6 


Bd ^ pO J/t 




4.44 


4.04-4.63 




5.05 


3.61 




1.43 


< 25 


B+ p+J/^ 




9.22 


8.39-9.62 




10.48 


7.49 




2.97 


< 77 


Bd -> K"K" 


0.09 


[0.01] 


0.05-0.56 [0.00-0.41] 


0.09 [0.005] 


0.08 [0.01] 


0.13 


[0.03] 


< 1.7 


Bd -> K+K' 


0.00 


[0.00] 


0.00-0.50 [0.00-0.50] 


0.00 [0.00] 


0.00 [0.00] 


0.00 


[0.00] 


< 1.7 


B+ K+RO 


0.09 


[0.01] 


0.05-0.58 [0.00-0.43] 


0.09 [0.005] 


0.08 [0.01] 


0.13 


[0.03] 


< 2.1 



Table 1: Predictions for yet unmeasured BRs. The channels above the horizontal line 
are those enhanced by charming penguins. In square brackets, the BRs obtained without 
charming penguins are shown. BR{B^ — ^ cuK^), BR{B^ ^t^^), BR{B^ rj'K^) 
and BR{Bfi r]'K^), for which preliminary measurements exist ('(LStgig ± 0.3) x 10~^, 
(1.1^^:^ ± 0.2) X 10-^ (7.ll|f ± 0.9) X 10-5 and (5.3l|^ ± 1.2) x 10"^ respectively), have 
been given in this table since they have not been used in the fit of tjl and 5l- 
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Channel 


QCDSR-CV 
BR X 10^ 


Lattice-CV 
BR X 10^ 


QM-CV 
BR X 10^ 


ABLOPR-CV 
BR X 10^ 


Exp. upper bound 
BR X 10^ 


Bd - vrO^o 


17 


11 


17 


4 


< 33 


Bd TT^D*^ 


28 


19 


35 


6 


< 55 


Bd - 


8 


10 


10 


3 


< 55 


Bd p^D*^ 


32 


40 


28 


19 


< 117 


Bd ^ ujD^ 


8 


10 


9 


3 


< 57 


Bd ^ u:D*^ 


32 


40 


27 


19 


< 120 


Bd ri'D^ 


8 


5 


8 


2 


< 33 


Bd ^ r?'5*o 


12 


8 


14 


48 


< 50 



Table 2: Predictions for emission- dominated decay channels for which only experimental 
upper hounds exist. 



Channel 


QCDSR-CV 


QCDSR-BV 


Lattice-CV 


QM-CV 


ABLOPR-CV 


Experiment 




BR X 10^ 


BR X 10^ 


BR X 10^ 


BR X 10^ 


BR X 10^ 


BR X 10^ 




301 


292-344 


293 


308 


318 


310 ± 44 


Bd TT + D*- 


323 


313-369 


314 


317 


341 


280 ± 41 


Bd ^ p+D- 


794 


770-907 


771 


802 


839 


840 ± 175 


Bd ^ p+D*- 


994 


964-1136 


966 


916 


1051 


730 ± 153 


B+ ^ TT+SO 


508 


498-527 


491 


498 


437 


500 ± 54 


B+ ^ n+D*o 


605 


595-619 


595 


604 


504 


520 ± 82 


B+ p+D» 


1015 


988-1103 


1078 


1022 


1015 


1370 ± 187 


B+ p+D*" 


1396 


1364—1496 


1553 


1269 


1479 


1510 ± 301 


Bd ROJ/^ 


81 


80-81 


65 


96 


103 


85± 14 


Bd -> K*°J/<ii 


164 


163-165 


185 


120 


51 


132 ± 24 


B+ K+J/^ 


84 


84 


68 


100 


107 


102 ± 11 


B+ K'+J/-^ 


171 


170-172 


192 


126 


54 


141 ± 33 


Bd ^ D+D- 


816 [880] 


710-1202 [769-940] 


801 [855] 


792 [864] 


935 [931] 


740 ± 284 


Bd ^ Dt+D- 


824 [887] 


723-1253 [844-996] 


807 [861] 


802 [871] 


842 [938] 


1140 ± 505 


Bd ^ D+D*- 


643 [691] 


565-994 [677-796] 


629 [671] 


626 [579] 


734 [731] 


940 ± 332 


Bd ^ Dt+D*- 


2358 [2535] 


2067-3584 [2412-2848] 


2309 [2462] 


1992 [2165] 


2693 [2682] 


2000 ± 729 


B+ ^ D+DO 


847 [913] 


736-1247 [797-975] 


831 [887] 


822 [896] 


970 [966] 


1360 ± 433 


B+ Dt+D° 


857 [922] 


752-1304 [877-1035] 


840 [895] 


833 [906] 


979 [975] 


940 ± 386 


B+ -> D+D*° 


669 [718] 


588-1034 [704-828] 


655 [697] 


651 [706] 


763 [760] 


1180 ± 462 


B+ D*+D*" 


2448 [2632] 


2146-3721 [2505-2956] 


2397 [2556] 


2067 [2247] 


2796 [2785] 


2700 ± 1045 


Results of the fit 




0.47 


(0.40-0.57) 


0.50 


0.50 


0.29 






0.42 


(0.21-0.68) 


0.00 


0.53 


0.00 




X^/dof 


1.00 


0.99-1.63 


1.64 


0.67 


1.97 





Table 3: Theoretical predictions for several BRs obtained by using different models for the 
semileptonic heavy-heavy and heavy-light form factors. For each model, ^ and 5^ are fitted 
by minimizing the /dof computed using the BRs above the horizontal line only. The fitted 
values of these parameters are also shown, together with the corresponding x^/dof. For 
QCDSR-CV, we also give in parentheses the value of ^ and obtained by putting 5^ = 0. 
Charming penguin contributions, using til and Sl as determined from BR{B^ — > K^tt'^) and 
BR{Bd — > K^t:~), are included. In square brackets, the BRs obtained without charming 
penguins are shown. 
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cos 5 




BR(B+ K"tt+) 
BR X 10'"'' 


BR{Ba -> K+TT-) 
BR X 10^ 


BR{Ba n+n-) 
BR X 10" 


VL 


Sl 


Experiment 






^••5_o. 9-0.2 "--^ 


+U.5+U.i , 1 


< 1.5 






QCDSR-CV 
LO = 5 GeV 


-0.7 
0.38 

+0.7 
-0.7 
0.38 

+0.7 
-0.7 
0.38 

+0.7 


0.48 
0.48 
0.48 
0.36 
0.36 
0.36 
0.24 
0.24 
0.24 


1.49 [0.18] 
1.62 [0.19] 
1.68 [0.20] 

1.36 [0.19] 
1.49 [0.19] 
1.54 [0.20] 
1.23 [0.19] 

1.37 [0.19] 
1.42 [0.19] 


1.61 [0.48] 

1.60 [0.19] 
1.59 [0.11] 

1.62 [0.40] 

1.61 [0.19] 

1.61 [0.12] 

1.63 [0.33] 

1.62 [0.19] 
1.62 [0.14] 


2.09 [0.93] 
1.59 [1.10] 
1.49 [1.15] 
1.30 [0.50] 
0.99 [0.63] 
0.94 [0.67] 
0.69 [0.21] 
0.52 [0.30] 
0.51 [0.32] 


-0.26 
-0.29 
-0.30 
-0.25 
-0.28 
-0.29 
-0.24 
-0.27 
-0.27 


1.73 
0.38 
0.06 
1.71 
0.32 

0.004 
1.62 
0.18 

-0.10 


QCDSR-CV 
NLO = 5 GeV 


-0.7 
0.38 
+0.7 
-0.7 
0.38 
+0.7 
-0.7 
0.38 
+0.7 


0.48 
0.48 
0.48 
0.36 
0.36 
0.36 
0.24 
0.24 
0.24 


1.43 [0.53] 
1.75 [0.56] 
1.85 [0.57] 
1.30 [0.53] 
1.62 [0.56] 
1.71 [0.56] 
1.22 [0.54] 
1.50 [0.55] 
1.57 [0.56] 


1.62 [0.97] 

1.58 [0.49] 
1.57 [0.35] 

1.63 [0.85] 

1.60 [0.49] 

1.59 [0.39] 
1.63 [0.75] 

1.61 [0.51] 

1.60 [0.44] 


2.24 [0.96] 
1.71 [1.26] 
1.64 [1.34] 
1.30 [0.51] 
1.04 [0.73] 
1.02 [0.80] 
0.47 [0.20] 
0.53 [0.35] 
0.54 [0.40] 


-0.34 
-0.28 
-0.26 
-0.32 
-0.25 
-0.24 
-0.25 
-0.22 
-0.21 


1.43 
0.25 

-0.08 
1.25 
0.15 

-0.16 
0.65 
-0.4 

-0.31 


ABLOPR-CV 
LO M = 5 GeV 


-0.7 
0.38 
+0.7 
-0.7 
0.38 
+0.7 
-0.7 
0.38 
+0.7 


0.48 
0.48 
0.48 
0.36 
0.36 
0.36 
0.24 
0.24 
0.24 


2.30 [0.54] 
2.30 [0.57] 
2.30 [0.58] 
2.03 [0.55] 
2.19 [0.17] 
2.23 [0.57] 
1 7s fn ^Rl 

1.93 [0.56] 
1.97 [0.57] 


1.50 [1.27] 
1.50 [0.52] 

1.50 [0.30] 

1.55 [1.05] 
1.52 [0.49] 

1.51 [0.33] 
1 "^Q fn S71 

1.57 [0.49] 

1.56 [0.38] 


4.67 [2.50] 
3.55 [2.97] 
3.41 [3.11] 
2.81 [1.36] 
2.37 [1.71] 
2.27 [1.81] 
1 41 fn ^fi 

1.21 [0.80] 
1.17 [0.87] 


1.07 
-0.40 
-0.43 
1.01 
-0.33 
-0.38 
0.93 
0.28 
0.16 


-0.23 
0.43 
0.14 
-0.23 
0.82 
0.50 
-0.25 
-0.89 
-1.35 


ABLOPR-CV 
NLO = 5 GeV 


-0.7 
0.38 
+0.7 
-0.7 
0.38 
+0.7 
-0.7 
0.38 
+0.7 


0.48 
0.48 
0.48 
0.36 
0.36 
0.36 
0.24 
0.24 
0.24 


2.29 [1.42] 

2.30 [1.50] 
2.30 [1.52] 
2.01 [1.43] 
2.30 [1.49] 
2.30 [1.51] 
1.76 [1.44] 
2.03 [1.48] 
2.09 [1.49] 


1.50 [2.50] 
1.50 [1.26] 
1.50 [0.89] 
1.55 [2.20] 
1.50 [1.27] 
1.50 [0.99] 
1.59 [1.93] 
1.55 [1.31] 
1.54 [1.12] 


5.09 [2.51] 
3.71 [3.29] 
3.60 [3.52] 
3.03 [1.34] 
2.56 [1.92] 
2.33 [2.09] 
1.49 [0.54] 
1.29 [0.92] 
1.25 [1.04] 


1.43 
-0.24 
-0.18 
1.34 
-0.32 
-0.22 
1.18 
0.40 
0.23 


-0.22 

0.64 
-0.06 
-0.23 
1.16 
0.51 
-0.26 
-0.68 
-0.91 



Table 4: BRs for several values of the CKM parameters cos 5 and a and different choices 
of the Wilson coefficients. The results are given for QCDSR and ABLOPR only. Results 
with LQCD and QM are very similar to those obtained with QCDSR. In square brackets we 
give the results obtained without charming penguins. 
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